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Abstract 

Microcavity polaritons are composite half-light half-matter quasi-particles, which have recently 
been demonstrated to exhibit rich physical properties, such as non-equilibrium Bose-Einstein con- 
densation, parametric scattering and superfluidity. At the same time, polaritons have some impor- 

Q 

tant advantages over photons for information processing applications, since their excitonic com- 
ponent leads to weaker diffraction and stronger inter-particle interactions, implying, respectively, 

O 

i tighter localization and lower powers for nonlinear functionality. Here we present the first experi- 

mental observations of bright polariton solitons in a strongly coupled semiconductor microcavity. 

> 

The polariton solitons are shown to be non-diffracting high density wavepackets, that are strongly 
localised in real space with a corresponding broad spectrum in momentum space. Unlike solitons 
known in other matter-wave systems such as Bose condensed ultracold atomic gases, they are 

O 

t-h non-equilibrium and rely on a balance between losses and external pumping. Microcavity polari- 

ton solitons are excited on picosecond timescales, and thus have significant benefits for ultrafast 



switching and transfer of information over their light only counterparts, semiconductor cavity lasers 

H 

(VCSELs), which have only nanosecond response time. 



1 



Non-spreading localised wave-packets or solitons occur in a variety of non-linear classical 
and quantum photonic and matter-wave situations, where they reflect the distinct physical 
properties of the system under consideration. Light only solitons were first discovered pQ 
and most actively researched in optical fibers [21 13]. Fibers provide a very good experimental 
realisation of the one dimensional (ID) Nonlinear Schrodinger (NLS) equation - a paradigm 
model in soliton physics, which has simple soliton solutions shaped by opposing dispersive 
spreading and nonlinearity induced self-phase modulation of optical pulses. Subsequently, 
the concept of optical solitons has been extended into the spatial domain, where diffraction 
of propagating beams can be compensated by nonlinear changes in the material refractive 
index (2J H] . Formally, anomalous group velocity dispersion is equivalent to diffraction: both 
require positive (self-focusing) Kerr nonlinearity to form either temporal or spatial bright 
solitons [2]. 

Recently, matter wave solitons have been demonstrated in Bose condensed atomic gases, 
where the dispersive spreading of the localized wavepackets induced by the kinetic energy 
of positive mass atoms is balanced by the attractive interatomic interaction [51 Ej. Inter- 
actions between the atoms used in condensation experiments are often repulsive. In order 
to create bright solitons in these cases, the effective atomic mass must be made negative, 
as demonstrated using optical lattices [7J. The same effect can be achieved for light in 
photonic crystals (8]. Strong exciton photon coupling in microcavities results in an unusual 
and advantageous lower polariton branch dispersion exhibiting regions of either positive or 
negative effective mass [HI [TO] depending on the values of transverse momentum P, [TU] . The 
transition from positive to negative mass is associated with the point of inflection of the 
energy-momentum diagram, shown in Fig. [T|a). The negative mass of polaritons coupled 
with repulsive polariton-polariton interactions favour the formation of the bright solitons 
which we realise in the present work. 

Recent discoveries made with microcavity polaritons include condensation [ITj H2] , vor- 
tices [131 E]) superfluidity [15] and follow a pathway similar to the one which led to the 
observation of coherent matter wave solitons [51 E] . The observation of low threshold bista- 
bility, polarization multistability [T6HT9] and parametric scattering [20J of polaritons have 
further prepared the necessary foundations for the realisation of half-light half-matter soli- 
tons. In contrast to atomic systems, polariton microcavities operate under non-equilibrium 
conditions. Atomic systems are typically described by conservative Hamiltonian models, 
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such as NLS or Gross-Pitaevskii equations, while the polariton system is intrinsically non- 
Hamiltonian. In the photonic context, this system is often referred to as dissipative, which 
implies the importance of losses, but also implicitly assumes the presence of an external 
energy supply [2TJ E2] • Theoretically, conservative half-light half- matter solitons have been 
predicted in Ref. [23], while dissipative bright polariton solitons in microcavities have been 
numerically studied in Refs. [2H [25] . 

OBSERVATION OF BRIGHT POLARITON SOLITONS 

Here we report bright solitons in a GaAs semiconductor microcavity (see Methods section 
for sample and experimental details). The point of inflection of the polariton dispersion is 
found at the in-plane polariton momentum k ~ 2fim~ 1 (see Fig. 1(a)), which corresponds 
to a group velocity of noninteracting polaritons de/dk ~ 1.8/im/ps, where e is the polari- 
ton energy. We conduct our experiments with polaritons having momenta above the point 
of inflection, where the polariton mass is negative, M = h 2 (d 2 e/dk 2 )^ 1 ~ —11.2 • 10~ 35 kg 
(—1.25 • 10~ 4 m e ) at k ~ 2.4/im -1 . As for the lossy and pumped Gross-Pitaevskii equations, 
see, e.g. Ref. [26] . reasonable estimates of the nonequilibrium soliton parameters can be ob- 
tained from the conservative limit of the system. We can estimate the polariton soliton width 
for a given polariton density, iV using the well-known expression for the healing length of a 
quantum fluid w = 2h(y/2MgN)~ 1 . This is obtained by equating the characteristic kinetic 
energy K of the dispersive spreading of the wavepacket with width w (K = h 2 /(2Mw 2 )) to 
the potential energy U of the polariton-polariton repulsion (U = gN, where g ~ 10 yueV/im 2 
is the two-body interaction coefficient). For the typical soliton potential energies realised 
in our experiment U ~ 0.3 meV (corresponding to N ~ 30/zm~ 2 ) from which we deduce 
w ~ 2/im. The short free polariton lifetime (~ 5 ps) means that solitons emerging from 
this balance will traverse distances of ~ 10/iin, before they dissipate. In order to sustain 
these solitons for longer we need to provide a continuous supply of energy. To sustain these 
solitons, energy is pumped into the microcavity using a CW pump beam, focused to a 70/im 
(FWHM) Gaussian spot (Fig. [TJb)). 

We highlight the important properties of dissipative solitons which justify the choice of our 
experimental parameters to realise soliton production. The polariton density varies strongly 
across the bright soliton profile reaching its maximum at the center, while the pump beam 
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profile tends to hold the density at a quasi- const ant level across the much larger pump spot 
[2UEH]- Furthermore, due to localisation in real space the soliton profile in momentum space 
is broad. As a result, solitons can only be expected under conditions when the pump state is 
unstable with respect to spatially inhomogeneous perturbations at momenta different from 
that of the pump [24"1 125] . If the microcavity is driven by a pump beam slightly blueshifted 
with respect to the unperturbed LP branch, the pump polariton field exhibits bistability as a 
function of the pump beam power and angle (Fig. flic), Fig. The polariton bistability is 
usually accompanied by parametric (modulational) instability, which was extensively studied 
in microcavities [TBI ¥2T\ [25] . This instability is a particular case of polaritonic four- wave 
mixing. Bright solitons are excited on top of the stable background of the lower branch of the 
bistability loop (Fig. [jjc)) and can be qualitatively interpreted as locally excited islands of 
the modulationally unstable upper branch solution [24"l 125] . Broadband four- wave mixing of 
polaritons, expanding well beyond the momenta intervals with the parametric amplification, 
enables coherent scattering from the locally perturbed pump to the continuum of momenta 
forming the soliton (Fig. [l](a)). Simultaneously, the soliton formation requires the transverse 
momentum of the pump k p , i.e. the incident angle of the laser beam, to be such that 
the effective polariton mass is negative, ensuring self-focusing of the repulsively interacting 
polaritons [23]. This mechanism can lead to self-localisation only along the direction of the 
pump momentum (ID solitons) [21]. Simultaneous localisation in two dimensions involves 
more subtle physics and the corresponding solitons exist only in a narrow range of pump 
intensities [2"5] . 

Fig. [T|c) shows the bistable dependence of the polariton emission intensity collected from 
our device at nearly zero transverse momenta (direction normal to the cavity plane) as a 
function of the pump momentum (k p ) at the fixed energy of the pump 1.5363eV, which is 
~ 0.3meV above the unperturbed LP branch at k p ~ 2.37/im -1 . With increasing k p beyond 
the bistable interval the intracavity polariton field increases abruptly over the excitation 
pump spot. This transition is also accompanied by the strong parametric generation of 
polaritons into the state with k ~ [T6l ED]- According to the above discussion and the 
predictions of [2U EH], bright polariton solitons should exist within the bistability interval 
(Fig- Be)). 

We excite bright polariton solitons using a ps writing beam (wb) focused into a spot (Fig. 
[lib)) of diameter in the range 7/iin to 15/iin, which is small compared to the diameter of the 
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FIG. 1: (a) Dispersion (energy-momentum) diagram of the lower branch polaritons and schematic 
representation of the soliton spectrum and the excitation scheme, (b) Schematic of soliton ex- 
citation in the microcavity structure. The CW pump with incident in-plane momentum parallel 
to the X direction is focused into a large spot. The pulsed writing beam, also incident along X, 
focused into a small spot triggers soliton formation, (c) Bistable polariton density as a function 
of the pump momentum; The arrow indicates the state of the system created by the pump before 
incidence of the writing beam, (d-g) Streak camera measurements of the soliton trajectories along 
the X direction excited under different conditions, (d) and (e) show the components of the soliton 
in TE and TM polarizations, respectively. This soliton is excited with a 7/xm wide writing beam 
with in-plane momentum the same as that of the pump, (f) is the same as (d), but the writing 
beam has half the momentum of the pump, (g) is the same as (f), but the writing beam width is 
14/um. 

pump beam of 70/mi. The writing beam is TE and the pump beam is TM polarized. The 
experimental arrangement of the pump and writing beams is shown schematically in Fig. 
1(b), with both beams incident along the X-direction. The unperturbed polariton density 
(before the application of the writing beam) was in the state corresponding to the lower 
branch of the bistability loop (see arrow in Fig. []Jc)) with no indication of the parametric 
generation of polaritons with momenta different from the pump. The bistability domain was 
scanned by changing the pump momentum, where it was found that the optimal conditions 
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for clear soliton observations exist close to the right boundary of the bistability interval 
shown in Fig. [l](c). The solitonic emission was collected along a 2/im stripe of the streak 
camera image (see Methods) parallel to the direction of incidence of the pump. The light 
was collected inside the finite interval of angles corresponding to momenta < k < k p , 
thereby avoiding collection of the reflected pump beam, which otherwise leads to detector 
saturation. Figs. [T|(d) and (e) show spatio-temporal traces of the intensities in TE and TM 
polarisations, for typical non-diffracting and non-decaying propagating wavepackets excited 
with the writing beam. Here the 7^m writing beam arrives at position X = — 20/im at time, 
t = where the transverse momentum is the same as that of the pump k w b — k„. 

We now present the experimental results, supported by theory in the next section, which 
prove that the wavepackets arise from bright polariton solitons propagating across the ex- 
citation spot. Firstly, we demonstrate that the velocity of the soliton is independent of k w b. 
The role of the writing beam is to create a local perturbation of the pump state, which in 
turn results in soliton formation due to scattering from the pump state (Fig. [jjb)). It has 
been shown numerically for linearly polarised polariton solitons that their velocity is close 
to the group velocity of the polaritons at the pump momentum [21] . To test this prediction, 
we changed the momentum of the writing beam to half of the pump momentum, k w b = k p /2 
(Figs. l(f,g). We observe a soliton velocity in Fig. [T]T) of ~ 1.68/im/ps, the same as that in 
Figs. [jjd) and (e), where k w t = k p . The independence of the polariton soliton velocity on 
k w b is in sharp contrast to that expected for conservative solitons, where the soliton velocity 
is solely determined by the momentum of the excitation pulse [2] . 

Secondly, we show that the size of the soliton is determined by the pump and cavity 
parameters and is independent of the size of the writing beam, as is true for other types of 
dissipative solitons [211 122]. The size of the excited wavepacket (w) in the soliton regime 
is expected to be fixed by the potential energy U of the solitons, where U is of the order 
of the pump energy detuning with respect to the energy of the unperturbed lower branch 
polaritons (w ~ 1/y/U as discussed earlier). This is illustrated in Figs. [i^a-f), which show 
the profiles of polariton wavepackets along their propagation direction (X) at different times 
and positions for writing beam sizes of ~ 7/im (Fig. |2^a-c)) and ~ 15/zm (Fig. [2]^d-f)). In 
the initial stage of the soliton excitation, these writing beams produce polariton wavepackets 
of very different widths (see Figs. [2^a) and (d)), which then quickly evolve into solitons of 
the same size ~ 5/im as shown in Figs. [21(b) and (e), as expected for soliton formation. 
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FIG. 2: (a)-(c) Experimentally measured spatial intensity profiles of a soliton created by the 7/um 
writing beam at different times, showing excitation and decay of the soliton; (d)-(f) Intensity 
profiles of a soliton created by the 15/im writing beam; (g,h) Dynamics of the peak intensity and 
width of the wavepackets excited by the 7[im writing beams with different powers as they propagate 
across the cavity. The writing powers where the above parameters remain quasi-constant in the 
interval of 20 to 40/Uin correspond to the formation of solitons. 

Thirdly, we demonstrate the unambiguous non-spreading and non-decaying features of 
the soliton wave packets. Figures [2^g) and (h) show the dependences of the intensity and 
width of the excited wavepackets, respectively, versus their position as a function of the 
writing beam power, P w f,. Within the range of P w b from 1.6 mW to 2.3 mW, the soliton 
intensity and width are nearly constant. Furthermore, the width versus position dependence 
(Fig. [2](h)) exhibits a plateau in the interval from ~ 20 to ~ 40/im. These signatures of 
polariton soliton formation are further discussed in the modeling section. At the edge of the 
pump spot, the pump intensity is insufficient to maintain the system within the bistability 
region and thus to sustain solitons. This is the main factor determining the practical extent 
of the soliton trajectory, and leads to spreading and dissipation of the excited wavepackets 
at times greater than 40-45 ps when the edge of the pump spot is approached, as shown in 
Figs. Ud-g) and Fig. gc),(f). 

We note that the intensity of the writing beam must be high enough to enable soliton 
switching. (Fig. |2](h)). As one reduces the pump intensity below the threshold of 1.3mW 
(Fig. |2](h)) for soliton switching, we observe an abrupt decrease in the intensity and an 
increase in width of the wavepacket. Further reducing the pump intensity P w b < ImW, only 
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FIG. 3: Two dimensional streak camera measurements of a soliton traveling across the microcavity 
plane. Experimental conditions are the same as in Figs. [2ja-c). Cross sections taken along the Y 
direction indicates a soliton full width at half maximum of 5/iin. 

quickly decaying non-solitonic wavepackets are excited. For pump momenta that are below 
the point of inflection, the writing beam triggers switching of the whole pump spot [I7J and 
no soliton formation is observed. 

Figure [3] shows two dimensional images of solitons. The soliton size along the Y coordinate 
perpendicular to the propagation direction (see Fig. [ijb)) at 30 — 40ps is ~ 5/iin, narrower 
than the initial beam size. This indicates suppression of diffraction and localization in the 
direction where polaritons have a positive effective mass. Whereas localisation in the X 
direction described above arises from the interplay between the negative effective mass and 
repulsive interactions, the origin of localisation along Y is different. It can be interpreted 
in terms of the phase dependent parametric nonlinearity and the interaction of propagating 
fronts, see Ref. [25] . 

We note that our observations have a very different physical origin to the triggered optical 
parametric oscillator observations of Ref. [2]. In that case, the system is already in a high 
density phase before a writing beam is applied [30] and furthermore the system was pumped 
at low k p in the region of positive effective mass, where soliton generation cannot occur. As 
a result, propagating excitations of the polaritonic condensate were studied, as opposed to 
the bright solitons reported by us (see Supplementary Information for further discussion). 
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NUMERICAL MODELING OF SOLITON FORMATION 



In order to provide further insight into our experimental observations, we have performed 
a series of numerical simulations. This is important since the previous theoretical studies of 
polariton solitons (2U [25] are restricted to the linearly polarized case, while in the present 
setting, with orthogonally polarised pump and writing beams, coupling between the two 
polarisations is important and leads to the formation of elliptically polarised solitons. The 
model we have used includes equations for the TE (polarisation parallel to the cavity plane) 
and TM polarised optical modes and for the respective excitonic fields (see Methods for 
details). 

Transforming from the laboratory frame to the frame moving with an unknown veloc- 
ity, which is consistently determined together with the soliton profile |24j . we have found 
branches of elliptically polarised solitons, see Figs. |3](a,b). These solutions retain their spa- 
tial profile and all other characteristics for indefinitely long times. The soliton branches with 
positive slope in Fig. E[a) represent propagating solutions, which are stable with respect to 
perturbations. The negatively sloped branches correspond to unstable solitons. These ideal 
solitons have been calculated on top of the spatially homogeneous pump and move with 
velocities ~ 1.7/im/ps, in good agreement with experiment. 

Another series of numerical simulations has been conducted to demonstrate the excitation 
of solitons with lps pulses on top of the gaussian pump beam, see Figs. |i|c,d). Fig. |3|c) 
(cf. Fig. [2^g)) and Fig. [4](d) (cf. Fig|2^h)) show changes of the peak intensity and of the 
width of wavepackets traveling across the microcavity for several intensities of the writing 
beam. The saturation of both parameters with increasing intensity of the writing beam 
(Fig. |4^c)) as well as the quasi-constant soliton width in the region from 5 /im to 30 fim of 
the propagation length (Fig. |4](cl) ) demonstrate the transition to the soliton regime, and is 
fully consistent with the experimental observations in Figs. [2|g-h). The observed saturation 
behaviour indicates that as for other dissipative solitons [2TJ E2] the soliton width in our case 
is fixed by the pump and cavity parameters[2H [25] and not by the writing beam powers. 
This is in strong contrast to conservative solitons, which for the given system parameters can 
exist with an arbitrary width determined by the density N initially induced by the writing 
beam. 
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FIG. 4: Numerically computed ideal solitons (a,b) and dynamics of the soliton excitation (c,d). 
(a) The full green line shows the bistable spatially homogeneous polariton state. The full black 
line and the dashed red line show the amplitudes of the TE and TM components of the elliptically 
polarized soliton, respectively, (b) Spatial profiles of the TE (full black) and TM (dashed red) 
soliton components. The inset shows the soliton profile when the light with k > k p has been 
filtered out. (c) and (d) Dynamics of the soliton width and peak intensity as functions of position 
across the microcavity. The writing beam is 7/iin wide and its duration is lps. Different line colors 
correspond to the different intensities of the writing beam indicated in the inset in (d). One unit 
of all the plotted dimensionless field amplitudes corresponds to the pump beam amplitude E p at 
the right end of the lower branch of the bistability loop, see (a). 

SOLITON OBSERVATION IN ENERGY-MOMENTUM SPACE 



The soliton formation can also be viewed as scattering of pump polaritons into a con- 
tinuum spanning a broad range of transverse momenta from to 2k p (see Fig. [TJa)) [24]. 
Our measurements of the energy-momentum (E(k)) profile of the solitonic wavepackets over 
time, as shown in Fig. [5j confirm this point. Initially, for t < 20ps, the polariton emission 
is concentrated mostly close to k p , whereas at later times (30 — 40 ps) the emission is dis- 
tributed over a broad range of momenta a characteristic feature of soliton formation [24J . 
The soliton spectra are expected to form a straight line tangential to the dispersion of non- 
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FIG. 5: Images of polariton emission in energy-momentum space taken at different times after 
the application of the writing beam in TM (bottom panel) and TE (top panel) polarisations, (a- 
j) Experimental data. Dashed line is a linear fit to the experimental measurements. The energy 
resolution of 0.1 meV limits the time resolution to ~ 10 ps. (k-n) Numerical simulations of polariton 
emission in time intervals of 13 — 27 ps (k,m) and 27 — 40 ps (l,n) after the application of the writing 
beam. Note, the span in the momentum space is wider than in the experimental figure. Dashed 
lines show the dispersion of the noninteracting LP branch. 

solitonic radiation [3JGI]. Indeed, the measured dispersion at 30 — 40ps can be approximated 
by a straight line (dashed lines in Fig. |5](b,c)) in agreement with our numerical modeling, 
see Fig. [H] (k,m). We note however that soliton emission forms at fc-vectors k > 0.5 /an" 1 
and spans over the point of inflection, where the dispersion of non-interacting polaritons is 
also close to linear. Nevertheless, the difference between the solitonic dispersion observed 
in experiment at around 30 — 40ps (in modeling at 13 — 27 ps) and the dispersion of the 
nonsolitonic polaritons observed at later times > 50ps (in modeling at 27 — 40ps) is clear 
from the data shown in Fig. JEJ The soliton velocity, deduced from the dashed lines in 
Figs. [5]^b, c) is ~ 1.6/mi/ps, consistent with the measurements in Fig. [TJand the numerical 
modeling. The difference between the experimental TE and TM emission in momentum 
space may arise from anisotropy of polariton scattering dynamics (see for example [3T]). 
which is not taken into account in our modeling. 
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DISCUSSION AND CONCLUSIONS 



The experimental observations of bright polariton solitons reported above open opportu- 
nities for the exploration of their potential applications in ultrafast information processing, 
since their picosecond response time is three orders of magnitude faster than that observed 
for the pure light cavity solitons in VCSELs j22j E2J [33] . Furthermore polaritonic non- 
linearities are 2-3 orders of magnitude larger than nonlinearities in VCSELs [22l 132] 133] . 
The measured transverse dimensions of polaritonic solitons are ~ 5/im (resolution limited), 
whereas the numerical model using the experimental parameters predicts ~ 2/zm, several 
times less than the 10/xm width typical for VCSEL solitons. Polariton solitons can poten- 
tially be used as more natural information bits than the propagating domain walls in the 
recently proposed integrated polariton optical circuits and gates (polariton neurons) [34]. 
The number of polaritons in the solitons observed is of the order of hundreds, which puts 
them into the category of mesoscopic structures. Realisation of spatially modulated ID 
microcavity structures, with reduced numbers of particles may create conditions for obser- 
vation of quantum solitons, and nonlinear functionality relying on few polariton quanta, 
potentially allowing designs of quantum polariton-soliton circuits. We note finally that un- 
like the well-studied polariton condensates which correspond to macroscopic occupation of 
a single state in momentum space, the highly occupied polariton soliton is strongly localised 
in real space with a broad spread in energy and momentum. 

In parallel with our work, Amo et. al. have very recently reported observations of dark 
ID polariton solitons [35]. These solitons are formed by polaritons with positive mass in 
collision with obstacles. A most important physical difference with the bright polariton 
solitons reported here arises from the fact that in the experiment of Ref. [35], the solitons 
are observed outside the pump spot, necessary to avoid locking of the dark soliton phase to 
the pump phase in their one-beam experiments. Such an arrangement selects conservative 
solitons propagating in the presence of the unbalanced absorption. 

METHODS 

We use a GaAs-based device with six 15nm thick quantum wells (QWs) grown by molec- 
ular beam epitaxy. The experiments were performed at a temperature of 5K. Propagation 
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of the wave packets in real time was recorded using a Hamamatsu Streak Camera with a 
time resolution of 2ps. An aspheric lens with a focal length of 4 cm and numerical aperture 
NA = 0.44, provided optical resolution of about 4-5 /im for both soliton excitation and 
imaging. Two-dimensional images of solitons at different times were reconstructed from an 
array of ID images versus position X such as in Fig. 1 obtained at different positions Y on 
the pump spot. Solitons were driven by a TM polarised CW pump beam and triggered by a 
pulsed TE polarised writing beam with duration of 5 ps. The pump and writing beams were 
directed to the microcavity through the same optical path using a polarising beam splitter. 
The orthogonal polarisations of the beams were chosen to ensure transmission of maximum 
available powers from each of the lasers to the sample. 

Numerical studies were performed using mean-field equations describing the evolution 
of slowly varying amplitudes of the TM and TE cavity modes and of the corresponding 
excitonic fields i/>tm,te' 

d t E TM - i^- {d 2 x + d 2 y ) E TM + ( 7c - iS c ) E TM = in R ^ TM + E p (x, y)e lkpX , (1) 

d t E TE ~ i^- {dl + flj) E TE + ( 7c - iS c ) E TE = in R ^ TE - iE wb (x, y, t)e lK ^ x , (2) 

OtlpTM + (je ~ iS e ) i>TM + "J [(1 + r)\lf) T M\ 2 + 2\lf) TE \ 2 ] ^TM - ^(1 - r^M^TE = ^rE T 0) 
dtlpTE + he - iS e ) ^TE + ~ [(1 + r)\iji TE \ 2 + 2\^ TM \ 2 ] i> TE - ^(1 - r)^ TE ^\ M = 1^^.(4) 

Here m c = 0.27 • 10 _34 kg is the effective cavity photon mass, HQr = 4.9867meV is the 
Rabi splitting, hq c = hq e = 0.2 meV are the cavity photon and the exciton coherence decay 
rates, S e = —1.84 meV, 6 C = —2.34 meV, g > is the nonlinear parameter, which can 
be easily scaled away, r = —0.05 parameterizes the nonlinear interaction between the two 
modes [36] • E p (x,y) is the pump amplitude with the momentum k p , the corresponding the 
angle of incidence 8 = arcsin[KA p /(27r)] and E w b is the writing beam amplitude. 

For the case of the homogeneous pump E p (x,y) = const and E w f,(x,y,t) = the soliton 
solutions are sought in the form E TM)TE = A TM ^ TE (x — vt)e lkpX , iPtm,te = Qtm,te(x — 
vt)e lhpX . The soliton profiles and the unknown velocity v are found self-consistently us- 
ing Newton-Raphson iterations. In our simulations of the soliton excitation the system of 
Eqs. ([l])-((4]) has been solved directly using the split-step method. 
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SUPPLEMENTARY MATERIALS: PROPAGATING POLARITON 

WAVEPACKET IN TRIGGERED OPTICAL PARAMETRIC OSCILLATOR 
CONFIGURATION 

As was discussed in the main text, optical bistability of the pump field together with 
the negative polariton effective mass at k p are necessary conditions for formation of a stable 
solitonic wave. In this case the soliton results from switching in the bistable region (k p ~ 
2.38/xmr 1 ) from the lower to the upper state locally in a region of a few microns. By contrast, 
if the whole pump spot is switched on, i.e the system is in the upper state in the bistability 
region or the k-vector of the pump is above the upper bistability threshold then a soliton 
can not be formed. This is demonstrated in the following experiment. 

Firstly, we set the system to the upper state (Fig. 1(c) of the main text), where a 
condensed 'signal' state is formed at k = 0, with long range spatial coherence over 30/zm, 
due to polariton-polariton parametric scattering from the pump state into 'signal' and 'idler' 
states at k S i gna i = and kidier = 2fc p . This corresponds to the polariton optical parametric 
oscillator (OPO), which was studied in Ref. pp. 

Secondly, we introduce a pulsed writing beam (wb) at k w b ~ 1.2/zm -1 focused to a small 
spot of ~ 7/im and observe propagation of the resulting wavepacket. The wavepacket is 
amplified due to polariton-polariton scattering from the switched on pump state and hence 
propagates macroscopic distances, superimposed on the signal condensate at k = 0. Such 
a process corresponds to the so-called triggered OPO (TOPO), firstly reported by Amo et. 
al. in Ref. [2]. 

The TOPO propagating wavepacket exhibits very different physical properties to the 
soliton. Fig. SI (a) shows the TOPO intensity as a function of time and position X (see 
Fig. 1(b) of the main text) for k p ~ 2.38/im -1 . It is seen that the TOPO wavepacket propa- 
gates over long distances of up to 50fim. However in contrast to the soliton propagation, it 
broadens significantly from 7^m up to 15 — 20/im within ~ AOps time and propagates at a 
speed 1.3 times slower than the observed soliton. Moreover the TOPO wavepacket intensity 
decays with time, again in marked contrast to soliton behaviour of Fig. 2(g). 

The TOPO wavepacket created by the writing beam can be considered as a condensate 
excitation, which has the properties of as a diffusive Goldstone mode [3J. There is competi- 
tion between polariton-polariton scattering from the pump to the signal state at k = and to 
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FIG. S 1: (a)- streak camera trace of propagating TOPO wavepacket as a function time and 
position recorded for k p ~ 2.38/im _1 and writing beam size 7/im. (b) - size (FWHM) of TOPO 
wavepacket recorded in panel ( Bbj clS cl function of position. 

the localised TOPO wavepacket at k w f,. Since the pump state is switched on over the whole 
excitation spot there are no excitations in the pump state itself. The TOPO wavepacket 
propagation speed is determined only by the writing beam k-vector k w j, ~ 1.2/im -1 and not 
by that of the pump as in the case of soliton. The TOPO wavepacket expands since the 
polariton dispersion at k w b ~ 1.2/im -1 is described by a positive effective mass. 

The time dependence of the TOPO wavepacket FWHM is expected to be strongly de- 
pendent on the initial size of the writing beam. It can be easily shown from solution of a 
time dependent Schrodinger equation that the FWHM of a noninteracting wavepacket with 
initial FWHM A with a dispersion described by positive effective mass m e ff is given by 
2./A 2 /4 + ^#^4t 2 . It is seen that the size of the wavepacket will increase with time more 
slowly for initially larger A. This is consistent with the above observation (wavepacket with 
initial size 7/xm is expected to spread up to ~ 20fim within 50 ps) and those in Ref. [2], 
where the TOPO wavepacket of a size 15-20 //m along the propagation direction is initially 
created and its size changes very little within 40-50 ps. Indeed, using the above formula we 
can estimate that it should increase only by 2 microns for m e ff ~ 10~ 34 kg. 
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